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ABSTRACT 

The 26 long-duration gamma-ray bursts (GRBs) with known redshifts form a distinct 
cosmological set, selected differently than other cosmological probes such as quasars 
and galaxies. Since the progenitors are now believed to be connected with active star- 
formation and since burst emission penetrates dust, one hope is that with a uniformly- 
selected sample, the large-scale redshift distribution of GRBs can help constrain the star- 
formation history of the Universe. However, we show that strong observational biases in 
ground-based redshift discovery hamper a clean determination of the large-scale GRB 
rate and hence the connection of GRBs to the star formation history. We then focus 
on the properties of the small-scale (clustering) distribution of GRB redshifts. When 
corrected for heliocentric motion relative to the local Hubble flow, the observed redshifts 
appear to show a propensity for clustering: eight of 26 GRBs occurred within a recession 
velocity difference of 1000 km s^^ of another GRB. That is, four pairs of GRBs occurred 
within 30 h^^ Myr in cosmic time, despite being causally separated on the sky. We 
investigate the significance of this clustering using a simulation that accounts for at least 
some of the strong observational and intrinsic biases in redshift discovery. Comparison 
of the numbers of close redshift pairs expected from the simulation with that observed 
shows no significant small-scale clustering excess in the present sample; however, the 
four close pairs occur only in about twenty percent of the simulated datasets (the precise 
significance of the clustering is dependent upon the modeled biases). We conclude with 
some impetuses and suggestions for future precise GRB redshift measurements. 

Subject headings: cosmology: miscellaneous — cosmology: observations — gamma rays: 
bursts 



1. Introduction 



A growing body of evidence suggests that long-duration GRBs arise from the death of massive 
stars. If true, the large-scale redshift distribution of GRBs should trace the star-formation history 
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of the Universe (Totani 1997; Wijers at al. 1998; Lamb & Reichart 2000; Porciani & Madau 2001; 
Bromm &: Loeb 2002). And, as explosive and transient events detectable to high redshift, GRBs 
clearly offer a probe of the Universe that is distinct from other lighthouses (Loeb 2002). Already, 
GRBs have helped shed light on the nature of damped Lyman a absorbers (e.g., Fynbo et al. 2002) 
and the nature of the faint end of the luminosity function of galaxies at moderate redshifts (e.g., 
Djorgovski et al. 2001b). 

As a redshift sample selected from an apparent isotropic population and from hosts spanning 
the breadth of the galaxy luminosity function, GRBs have the potential to provide a unique census 
of the redshift distribution of matter in the Universe. Not only can the large-scale distribution of 
GRBs help constrain the star-formation rate (SFR), but, in a manner complementary to pencil- 
beam surveys and magnitude-limited galaxy and quasar surveys, the small-scale distribution can 
help test hypotheses and observational suggestions about clustering and periodicity of sources in 
redshift space. 

Given the penetrative powers of GRBs through dust, several studies have attempted to compare 
the large-scale redshift distribution of GRBs with the SFR obtained by other means (Schaefer et al. 
2001; Lloyd-Ronning et al. 2002a; Stern et al. 2002; Norris 2002). Such studies focused on using 
the prompt 7-ray properties of BATSE bursts, calibrated with some measured GRB redshifts, to 
determine the bursting rate. With so few actual redshifts measured, each of which were measured 
under different observational constraints, the task of distinguishing different SFR scenarios from 
the sparsely sampled GRB rate is considerably challenging. Indeed we show herein, that when the 
observational biases in redshift determination are taken into account, several proposed forms for 
the universal SFR cannot be distinguished. 

Investigation of the small-scale distribution of GRBs has not been presented thus far. Of 
particular interest is whether the clustering properties in redshift space are significant in light 
of the observational biases and the intrinsic redshift distribution. Do GRBs occur at preferred 
redshifts? There is enough historical, albeit controversial, evidence to warrant an investigation 
of the significance of special redshifts in the Universe using GRBs. Broadhurst et al. (1990), for 
instance, found evidence for periodicity of clustering of redshifts on 128 hT^ Mpc scales in two 
pencil-beam surveys. While excess power on these scales has not been definitively confirmed in 
higher-redshift studies (e.g., with the 2dF redshift survey: Hawkins et al. 2002; although see Gal h 
Djorgovski 1997 and Duari et al. 1992) the Broadhurst et al. result may be statistically significant 
in comparison with A'^-body simulations (Yoshida et al. 2001). But, as Yoshida et al. emphasize, 
large clusters in a galaxy sample tend to accentuate the appearance of excess power at certain 
scales^. A GRB sample docs not suffer such a bias since bursts occur in regions of space that are 



^Excess power on such 128 hT^ Mpc scales could be an imprint of a primordial density fluctuation (Dekel et al. 
1992). As a statistical measure, if this interpretation is true, then the same redshifts with over-densities would not 
reoccur in surveys that sample different directions on the sky. However, specific special redshifts would be preferred 
if the observed power arises by more exotic means. For example, oscillations of the scalar potential (Morikawa 1991) 



-3- 



gravitationally (and causally) disconnected. 

The focus of this paper is on the small-scale (clustering) distribution of GRB redshifts. In 
addition we examine the large-scale distribution of GRB redshifts, and note the difficulty in de- 
termining the global GRB rate due to the small sample and strong observational biases. In §2 we 
present the GRB rcdshift sample for 26 bursts, corrected for the heliocentric motion through the 
local Hubble flow and construct an observationally-motivated probability distribution for redshift 
discovery. In §3, we compare the large-scale distribution of GRB redshifts with that expected from 
the observing biases and an intrinsic rate distribution. After showing the insensitivity of the present 
sample to distinguishing various SFRs, and under the anzatz that GRBs should trace the SFR, we 
fix a model for redshift discovery probability that adequately reproduces the universal SFRs. Using 
this distribution, we then test the significance of the number of observed pairs of GRBs with small 
recession velocity differences. We conclude with a discussion about the importance and usefulness 
of conducting precise redshift measurements for future GRBs. 



2. The Redshift Sample 

Table 1 presents a summary of the sky location and observed redshifts associated with 26 
cosmological GRBs, in order of increasing redshift. The highest observed redshift system derived 
using absorption or emission features with the highest reported accuracy is given in column seven. 
The line features used to measure the reported redshift are given in column six. Absorption-line 
redshifts, which place a strict lower limit to the actual GRB redshift, are found when a spectrum of 
the afterglow is absorbed by metal-line systems along the line of sight. Emission-line redshifts are 
derived from spectroscopy of the galaxy associated with the GRB position, either found promptly, 
superimposed with the spectrum of the transient afterglow (e.g., 980703; Djorgovski et al. 1998), 
or at later times, when the afterglow light has faded. Nineteen sources in the sample have mea- 
sured emission-line redshifts, and 12 have absorption-line redshifts. The uncertainty in the rcdshift 
measurements are usually provided in the literature, but where none were provided we assumed an 
error of 10 percent in the least significant digit reported. 

Since we are interested in testing the observed rate and redshift distribution of GRBs against 
different models for star-formation and small-scale clustering, we now discuss the observational 
biases in GRB detection, localization, and finally, redshift determination. 



or the gravitational constant (Salgado et al. 1996) induce an oscillation in the Hubble constant versus cosmic time. 
The result is an apparent clustering of sources in redshift space at epochs where the expansion temporarily slows. 
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2.1. Biases in GRB Detection 

If GRBs arise in coUimated emission, then the relativistic Doppler beaming ahows for detection 
of a GRB only if the detector lies in the direction of the cone of emission. By studying the opening 
angles of those GRBs which have been detected, it is reckoned that only about one in every five 
hundred GRBs are detectable at (i.e., pointed toward) Earth due to beaming (Frail et al. 2001). 

Of those bursts that are beamed toward Earth, only those that reach a certain critical flux level 
will be detected by the triggering algorithms of a given satellite. The typical threshold for on-board 
detection in the BATSE catalog is a peak flux of a 0.3 photon cm~^ (in 't Zand & Fenimore 1994; 
Pendleton et al. 1998). In finding that the total energy release in GRBs is approximately constant, 
Frail et al. (2001) showed that bursts that are more highly collimated appear to be intrinsically 
brighter (in flux/fluence per unit solid angle). Thus, for a scenario where GRBs are jetted and 
release about the same amount of electromagnetic energy, brighter bursts tend to be more distant. 
As a clear illustration of this, note that the peak flux of the highest rcdshift burst, GRB 000131 
{z = 4.5), was at the top 5% of the entire BATSE catalogue (Andersen et al. 2000). 

This somewhat counterintuitive trend implies that the brightest GRBs can be detected to 
extremely high rcdshifts, certainly beyond the age of reionization {z ~ 7—10). More important for 
the present study is that Malmquist bias should play little role in diminishing the observed burst 
rate at high-redshift: bursts significantly fainter (by up to ~ 10~^) than GRB 000131 could have 
been detected at similar redshifts. 

Using the prompt high-energy properties, irrespective of individual burst redshifts, it has been 
shown that the brightness distribution of GRBs is consistent with a broad range of star-formation 
rates and luminosity functions (e.g., Loredo &; Wasserman 1998; Kommers et al. 2000; Schmidt 
2001; Stern et al. 2002). Using distance indicators calibrated from a sub-set of the bursts with 
known rcdshift, Lloyd-Ronning et al. (2002b) claimed that the GRB rate increases monotonically 
at least out to z ~ 10. In our opinion, this suggestion points more to the pitfalls of extrapolating 
GRB distance indicators beyond the redshift range in which they were calibrated, than revealing 
concrete properties of the high-redshift bursting rate density. Another approach to constrain the 
burst rate, which we attempt herein, is to use only those bursts with actual measured redshifts. 

2.2. Biases in Afterglow Localization 

For those GRBs that are triggered and rapidly localized using the prompt X-rays or 7-rays, the 
observing community tries to determine a sub-arcsec position and ultimately associate a redshift via 
host or absorption-line spectroscopy. This is accomplished by first identifying a transient afterglow. 

Three major biases are important for afterglow discovery. First, GRBs will be preferentially 

localized if they occur at a time and place in the sky where ground- or space-based telescopes can 
rapidly observe the prompt burst position with large enough fields-of-view. Second, only GRB 
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afterglows with brightnesses above a given threshold for the particular afterglow-discovery observa- 
tions will be localized. Last, only GRB afterglows that are not heavily extincted by dust obscuration 
will be first localized at optical wavelengths. Clearly, radio and X-ray afterglow discoveries do not 
suffer this third bias. The existence of a population of dark bursts (Djorgovski et al. 2001a; Piro 
et al. 2002), those bursts without detectable optical afterglow emission, is clearly a result of all 
three biases. (The current dark burst debate resolves around the relative importance of dust versus 
observational biases in the manifestation of dark bursts (Berger et al. 2002; Fynbo et al. 2001).) 

Since, for now, GRBs with redshifts are preferentially those that were first localized at optical 
wavelengths (all but 970828, 980329, and 000210 were first localized at optical wavelengths), it seems 
appropriate to suppose that a population of optically-localized GRBs should trace a star-formation 
rate derived from other optically-selected samples; Porciani & Madau (2001) have provided a pa- 
rameterized version of a SFR based upon the un-dust corrected Hubble Deep Field measurements 
(model SFl); in this model the SFR drops beyond 2; ^ 2. In future GRB redshift samples of bursts 
that are first localized at radio. X-ray or even infrared wavelengths, the expectation is that the 
observed rate will more closely trace the true (high mass) SFR in the universe. 

2.3. Biases in Redshifts Determination 

For those GRBs that are triggered and localized, most are followed up spectroscopically. If 
an absorption or emission-line redshift is found, it is of interest to know the relationship between 
this and the true GRB redshift. A definitive measure of burst redshifts would need to connect the 
bursting location with some stationary gas local to the GRB, resulting in a detection of transient 
features. To date, no transient absorption or emission features at optical wavelengths features have 
been associated with a GRB afterglow. Transient features, consistent with (later-time) optical 
spectroscopic redshift measurements, have been seen in a prompt burst spectrum (Amati ct al. 
2000; Le Floc'h et al. 2002) and a handful of X-ray afterglow spectra (e.g., Piro et al. 1999, 2000). 
At present, the uncertainties in these X-ray-derived redshifts arc orders of magnitude larger than 
those found typically at optical wavelengths. Perhaps more important, redshifts derived purely from 
X-ray spectroscopy are contaminated from an unknown (potentially relativistic) outflow speed of 
the emitting or absorbing material; thus, precise measurements of the systematic redshift of GRB 
progenitors are best determined using the UV and optical features of host HII regions and/or 
intervening host clouds. 

If a GRB occurs in intergalactic space, with no afterglow absorption due to local gas, then 
the observed redshift would be systematically lower than the true redshift. However, to date, the 
nature of the highest z absorption lines, suggests the presence of gas columns and abundances not 
typically associated with tenuous outflow gas halos of galaxies. Specifically, the equivalent widths of 
the highest redshift absorption systems in GRB afterglows are significantly larger than in systems 
seen though quasar lines-of-sight (e.g., Salamanca et al. 2002). This is naturally understood in the 
context of the prevailing progenitor model, namely that long-duration bursts arise from the death 
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of massive stars and, owing to the rapid evolution of massive stars, GRB explosion locations should 

be in or near relatively dense regions of active star-formation. Thus, GRBs should reside in or near 
the region that gives rise to the gas absorption. Consistent with this picture is the observation that 
well-localized GRBs, a subset of which comprises the present sample, have been observationally 
connected to the location of the light of the putative host galaxies (Bloom et al. 2002). 

The probability of a spurious (spatial) association of an afterglow with its putative host has 
been shown to be small (^ 10"^) for most bursts. In fact in a sample of 20 GRBs with arcsecond 
localizations, statistically, at most a few spatial associations could have been spurious (Bloom 
et al. 2002). Therefore, from the location argument, GRB redshifts determined from emission line 
spectroscopy of the hosts are likely to be close to the true systematic velocities of the progenitor; 
at most, we would expect a velocity offset on the order of a few hundred km s~^ from the motion of 
the progenitor system about its host. Likewise, unless the afterglow light pierces some high-velocity 
outflowing material, absorption-line redshifts should also provide a accurate measure of the true 
GRB redshift. Perhaps most convincing that measured redshifts are closely related to the real 
systemic redshifts, is that of the five GRBs with both emission and absorption line spectroscopy, 
the highest redshift absorption system is always found to be consistent with the emission-line 
redshift of the host. 

2.4. Correcting the measured redshifts for heliocentric motion relative to the local 

Hubble frame 

One striking feature of the current redshift sample is the small apparent redshift differences, on 
the order of a thousand km s~^, between some bursts from different sky locations and apparently 
random trigger dates (differences of 6 months to 5 years). (Unlike in pencil-beam surveys the 
proximity in redshift can have nothing to do with systems, such as galaxies in clusters, that are in 
causal contact.) These small velocity differences are comparable to the velocity of the Solar System 
through the local Hubble frame, and so a correction for this systematic motion is warranted. Such 
motion was measured by the angular dependence of temperature variations in the Cosmic Microwave 
Background (CMB) with the COBE satellite (Smoot et al. 1991; Fixsen et al. 1994). This measured 
velocity {Vq = 365 it 18 km s~^) and direction (toward Galactic coordinates / = 264.4 ± 0.3 deg, 
b = 48.4 lb 0.5 deg; Fixsen et al. 1994) are due to vector sum of the heliocentric motion about the 
Galactic center, the peculiar motion of the Milky Way in the Local Group and the Local Group 
in- fall. We assume that all reported redshifts have been corrected to the heliocentric redshifts; this 
correction is at most of order tens of km and so this assumption is relatively unimportant. 

In order to remove this systematic bias from a measured redshift z, we find the angle (Ocmb) 
between the GRB position and the heliocentric motion through the CMB; we then find the projected 
heliocentric velocity toward the GRB as f0,proj = VqCOS^cmb- If uc is the uncorrected apparent 
recession velocity of a distant source, then, from the definition of redshift and the formula for 
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relativistic velocity addition, the corrected redshift in the local Hubble frame (Ihf) is found from: 

1-,,^ n2 1 + TO0,proj/c + U + t;0,proj/c 
1 + Ui;o,proj/c -U- Vo^proj/C 

with the uncorrected apparent recession velocity as, 

u-ii±^ (2) 

The calculated values of z\^i are given in column eight of Table 1. The associated uncertainties are 
derived by error analysis of equation 1, assuming that the errors on Vq and ^cmb are uncorrelated 
and noting that the fractional error on ^cmb is significantly smaller than the fractional error on 
As seen in Table 1, for all but the most accurately measured redshifts, this correction adds 
negligibly to the fractional error of the redshift. 



3. Results 

Table 2 shows the absolute difference in redshift and recession velocity of each burst matched 
with the burst that is closest in apparent recession velocity. Figure 1 shows this distribution of 
"nearest neighbor" apparent recession velocity differences versus redshift. The smallest groupings 
occur around redshift of unity, where the density of known GRB redshifts is highest. Four pairs 
are closer than ~1000 km s~^ (0.0034 c) in redshift space^, and the largest velocity difference is 
between the two highest redshift bursts {\^v\ = 0.22c). 

Interestingly, the velocity differences of the four closest burst pairs (at redshifts z =0.692, 
0.844, 0.961, and 2.035) are less than the velocity dispersion of a large cluster of galaxies — and 
so, given the unknown peculiar velocity of the GRB hosts and the progenitor system within the 
hosts, these close burst pairs are consistent with having occurred at the same cosmic time. Even 
assuming no peculiar velocities of GRB progenitors relative to their local Hubble flow, eight of 26 
GRBs occurred within 30 Hq^ Myr of another GRB. This was determined by computing the look- 
back times between the four close pairs and assuming a cosmology with Hq = 65 km s~^ Mpc~^, 
Qa = 0.7 and J^m = 0.3. This cosmology is assumed throughout the paper where needed. 



3.1. Simplistic probability calculation 

How significant is this proximity in cosmic time? The bursts listed in Table 1 have been de- 
tected with recession velocities in the range wiow = 0.30 c (GRB 011121) to t^high = 0.94 c (GRB 000131) 



^Of the eight bursts in the close four pairs, six redshifts were determined via emission hues from the putative host 
galaxies. One redshift (GRB 000926) was one with absorption spectroscopy only and another redshift (GRB000301C) 
was found with both emission and absorption spectroscopy. 
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In an ensemble of n bursts, assuming uniform probability of a GRB occurring and being detected 

between vi^^ and t^higln the probability that given pair of successive bursts are not as close as a 
distance Avc is P = exp (— n (fhigh — ^^iow)/^'i'c)- Therefore, the probability that at least one close 
pair exists in the ensemble, is, 

P(l) = 1 - P"-i (3) 

For A^c = 1000 km s~^ and n = 26, the probability of at least one close pair with Av < 1000 km 
is 0.967. The probability of at least one close pair with at most the smallest velocity difference 
in the observed GRB sample {Avc = 259 km s~^) is 0.585. Both these probabilities thus indicate 
that there is nothing particularly unusual about the occurrence of at least one close pair in the 
current sample. 

In general, the probability that exactly k pairs are close can be approximated as, 

^ [ Av, \ {n-k)\kV ^ ' 

For the present sample, there are A; = 4 close pairs. Using this simplistic calculation we therefore 
expect four pairs to occur by random chance in 18.6% of samples. We have verified this approxima- 
tion by a Monte Carlo simulation, selecting n bursts uniformly over the range [fiow) I'high] i however, 
this approximation breaks down for large velocity differences as it does not account for relativistic 
velocity subtraction between successive bursts in velocity space. 



3.2. Probability calculation with a model distribution 

The above probability estimate does not take into account some important observational and 
endemic biases. First, the rate of GRBs is not uniform in redshift or recession velocity space. 
Instead, we expect GRBs to trace (or at least approximate) the star-formation rate (Totani 1997; 

Wijers et al. 1998), so that the peak of bursting activity should occur around z ~ 1 — 2. This will 
make close velocity pairs more likely at redshifts near unity. Second, the chance that a burst will 
be localized well enough to follow-up with spectrometers is not uniform in redshift. Instead, this 
chance depends, for optical localizations, sensitively on the native dust obscuration of the afterglow 
(see §2.2). Third, observing conditions and instrumentation play a strong role in determining 
whether the redshift of a GRB can be detected at that redshift. Most important, emission lines 
and absorption lines falling outside the broadband spectral coverage hamper the ability to detect 
redshifts (this may be an explanation for why no redshifts have been detected for 980519, 980326, 
and 980329). Even if a redshifted line falls within the spectral range, the presence of night sky lines 
make detection of redshifts at certain wavelengths more unlikely. Fourth, instrumental sensitivity 
varies across the spectral range and from instrument to instrument, night to night, and airmass to 
airmass. 

Since we wish to know how the observed redshift sample compares with that expected given 
the above considerations we construct a toy probability model, P{z)^ giving the differential prob- 
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ability that a burst at redshift z could ultimately yield a redshift measurement. We construct a 
global (rather than individual) distribution, neglecting the instrumental sensitivity differences from 
burst to burst and the non-negligible chance that a redshift will be incorrectly assigned even after 
the detection of one or more spectral features. We use several parametrized versions of the star 
formation rate (SFR) from Porciani h Madau (2001). The overall shape of the distribution, Pg{z)i 
is then proportional to the number of bursts per unit redshift per unit observer time; that is, Pg[.z) 
is proportional to the co-moving volume element and the co-moving SFR (psfr) versus redshift: 

D . s diV psFR dV 

^ ( exp(3.4.) \ ^ dl{z) 

Vexp(3.8z)+45; (1 + zf^^^{\ + zf + Ok(l + + ' 

following eqs. [2]-[4] in Porciani & Madau (2001). Here d^iz) is the luminosity distance and 
r^in + + = l- The factor of (1 + z) in the denominator accounts for the time dilation of 
the co-moving GRB rate. Following the discussion in §2.2, we have nominally chosen the un-dust 
corrected SFR from Porciani h Madau (2001) (SFl). We also explore the two other models for 
Pg(z) in Porciani h Madau (2001). In our model we do not take into account any biases related 
to the distance of the source in localization (in §2.2 we argued that there should be little bias 
in triggering based on distance). However, because of the observed anti-correlation between jet 
opening angles and prompt emission fluence (Frail et al. 2001), and the theoretical correlation 
between prompt burst emission and afterglow luminosity (Panaitescu et al. 2001), aside from dust 
obscuration, the afterglow from triggered GRBs that originate from higher redshift GRBs should 
not be significantly dimmer at the same observer time (see also Ciardi h Loeb 2000; Lamb h 
Reichart 2000). 

Nevertheless, if there is any dependence of triggering/localization on d^, we can try to account 
for such using a probability function related to luminosity distance, Pl{z)- There is also a clear 

effect of distance upon redshift discovery: once a burst afterglow is localized and a spectrum is 
acquired, the detectability of emission features is diminished with increasing luminosity distance 
(nominally as d^^). Systematically, only GRB hosts with higher rates of unobscured star formation 
will be detectable at higher redshifts for a given integration time. For a given afterglow brightness, 
the distance bias does not exist for absorption-redshift GRBs as long as the particular redshifted 
line is observable in the spectrum. To try to account for these elusive effects, we take the probability 
of redshift discovery due to distance [Pl(z)] as unity from z = to zi, decreasing with di to some 
power L. Nominally we take the values of zi = 1 and L = —2 but allow these quantity to vary in 
our modeling (see §3.3). 

Using the observability of emission and absorption lines in the spectral range and in the 
presence of night skylines, we construct a relative probability of redshift detection, Ps{z)- We take 
the spectrum range as 3800-9800 A and assume a redshifted line is not observable if it falls within 
a wavelength that is half the instrumental resolution of a strong sky line listed in Osterbrock & 
Martel (1992). Here, we assume the instrumental resolution of 5 A (dispersion of ~2. 5 A/pixel). 
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Based upon Table 1, we assume that a redshift can be obtained unambiguously with the detection 
of at least one of four star-formation emission features: Ly a A1216 A, [O II] AA3727 A, Ha 6563 
A, and [O III] AA 4959, 5008 A. Also based on Table 1, for redshifts based upon absorption features 
only, we require that at least three of the following lines are detectable: Fe II A2344.2 A, Fe II A 
2374.5 A, Fe II A 2382.8 A, Fe II A 2586.7 A, Fe II A 2600 A, Mg II AA 2796.4, 2803.5 A, Mg I A 
2853.0 A. If no redshifted star formation line nor prominent absorption line is observable at a given 
redshift, zq, then wc set Ps{zo) = 0. Even if a line is observable, it may be too weak in emission or 
too low in equivalent width to be detected. Therefore with less and less lines observable at a given 
z, the value of Ps{z) should be diminished. 

Figure 3 shows our toy model for Ps{z), based upon the above considerations. We stress that 
this is a toy model for ground-based optical spectroscopy. Different instruments will provide free- 
spectral ranges that differ from our nominal range (3800-9800 A) . The resolution of each instrument 
also will also particularly affect the observability of faint, narrow emission/absorption features in 
the presence of sky lines. A moderate-resolution spectrograph {R ^ 10000) should allow for feature 
detection closer in wavelength to a sky line; even features which partially overlap a skyline could 
be detectable (e.g., Bloom et al. 2003). 

The resultant relative probability distribution of detecting an optical redshift for a triggered 

GRB is constructed as P{z) = Pg(z) x Ps(z) x Pl{z) and is depicted in Figure 4. The relative 
probability distribution was calculated in bins of 6z = 0.001. As can be seen in the normalized 
cumulative distributions shown in Figure 2, the overall shape of the cumulative distributions are 
largely unaffected by the sky lines. The most pronounced feature in P{z) is the drop in detection 
probability from z 1.5 — 2, due to the inaccessibility of strong emission lines in the optical 
bandpass. The relative probability for detection of bursts at the redshift of GRB 980425 {z = 
0.0088) is exceedingly small and thus serves to emphasize the distinction between GRB 980425 and 
the other long-duration GRBs with known redshift (e.g., Galama et al. 1998; Bloom et al. 1998; 
Schmidt 2000; Porciani & Madau 2001). As such, we have not included GRB 980425 in our analysis. 

3.3. Consistency of the Large-Scale GRB Distribution with a Variety of SFR 

Models 

As shown in Figure 2, the large-scale shape of the observed redshift distribution is adequately 
described by the model for L = —2, with the KS probability that the observed deviations are 
consistent with a random selection from the model of Pks = 0-64 (SFl), 0.51 (SF2), and 0.32 
(SF3). Models where no account for the observational bias of detecting emission lines from high- 
rcdshift GRB hosts (L = 0) are clearly ruled out. While there are small differences in the KS 
probability between various models for the star-formation rate, such models cannot be statistically 
distinguished by the current sample of GRBs with known redshifts. However, given the low KS 
probabilities for L = 0, the data do support the notion that the observational biases of detecting 
emission lines from high-redshift GRB hosts must be taken in account. Note that a fairly large 
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range of L values (r^ — 1 to —3 for SFl) yields a KS consistency with the data. 

How sensitive are these results to our toy model for Ps{z) and Pziz)"^ Changing the value of 
zi we still get acceptable (but lower) KS statistics for L = -2: Pks = 0.20 {zi = 0.5; SFl), 0.22 
{zi = 1.25; SFl), 0.04 (z/ = 1.5; SFl). With zi > 4.5, the L = -2 case is effectively equivalent to 
the L = case {Pks = 0.003 for zi = 5; SFl). By removing the effect of the sky lines and limited 
free-spectral range in redshift determination, the KS statistics arc still acceptable, albeit with lower 
values of the KS statistic: Prs = 0.37 (SFl), 0.19 (SF2), 0.10 (SF3). 

3.4. Testing the Small-Scale GRB Redshift Distribution 

Assuming that GRBs do trace the global star-formation rate we fix the form of Pl{z) that 
gives a reasonable agreement with all three SF models. From here, we can test the significance 
of the small-scale clustering properties. To do so, we produce a Monte Carlo realization as sets 
of GRB redshifts drawn from the probability for redshift detection P{z) for each SFR model. We 
simulated 5000 iterations of sets of GRB redshifts. For each iteration, 26 bursts were selected 
uniformly from [0,1] and then mapped to redshift using the cumulative and normalized distribution 
of P{z), interpolating between bins to increase the resolution. 

We thought of no obvious existing statistic to compare the small-scale redshift structure of 
the observed distribution with the Monte Carlo set. However, since one feature is the existence of 
such close redshift pairs we can ask how often such numbers of pairs are found in the distribution. 
Table 3 summarizes the results of the comparison. Column two lists the number of observed GRBs 
within an apparent recession velocity of \^v\ of another GRB. Columns three, four and five give 
the probability of such occurrences in the simulated distributions for the three different GRB rates, 
Pg{z)- Following the table, the simulation predicts that at least two bursts (one pair) should occur 
by random chance in 53% of real world samples for the smallest observed velocity difference. This 
is close to the probability obtained in our simplistic calculation (59%) in §3.1. For Av values near 
1000 km s~^, the probability drops to f^l9% before rising again at large velocity offsets. Therefore, 
despite the apparent close pairings in redshift space, we find no significant small-scale redshift 
clustering in the present sample. 

Note, however, this comparison only references 2-point correlations. Following column eight of 
Tabic 1 and column 6 of Table 2, there arc two groupings of three bursts within 2500 km (at 
z = 0.69 and z = 0.84). Using the Monte Carlo set, the probability of getting 2 groupings of three 
bursts within 2500 km s"^ is 0.24 (SFl), 0.17 (SF2), and 0.13 (SF3). Again, the present sample 
shows no evidence for significant multi-burst clustering. 

With the advent of Swifi^ , it is not unreasonable to expect upwards of one hundred new GRB 



^http : //miv . swift . psu . edu/ 
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redshifts within the next several years. As the density of redshifts increases, the number of close 
redshift pairs (and triples) will certainly increase. Commensuratcly, for an expanded sample to yield 
significant clustering results, the number of close pairs must increase. Using the present models, as 
a function of bursts in a sample, we predict the required number of observed close pairs in order 
to be considered a significant excess above a random sample. Figure 5 depicts the results. For 
the present sample, approximately 1 (2) more close pairs would be needed with velocity differences 
less than 1000 km s~^, for the clustering results to be considered significant at the 0.05 (0.01) 
level. From the Figure, we see that in a sample of 50 GRB redshifts, there is a 5% random chance 
of getting 5 close pairs within about 250 km s~^. If 6 or 7 more pairs are found, this would be 
considered a significant result. 



4. Conclusions 

We have demonstrated, via a KS test, that the large-scale distribution of observed long-duration 
GRB redshifts is compatible with having been drawn from a reasonable model for the detection 
probability of burst redshifts. This supports the suggestion, based upon 7-ray properties, that the 

rate of GRBs rises rapidly out to redshifts of order unity (cf. Schmidt 2000; Lloyd-Ronning et al. 
2002a; Stern et al. 2002). To our knowledge, this is the first investigation to claim this using only 
the observed GRB redshifts and a model for the observational and intrinsic selection function for 
GRB redshift discovery. 

The data are relatively insensitive to various forms of the underlying star-formation for a 
given set of redshift selection functions Ps{z) and Pl{z). We have shown that by not considering 
the biases of ground-based optical spectroscopy in redshift determination (i.e., setting Ps{z) = 1) 
the large-scale redshift distributions are still consistent with the SFR expectations, although at a 
lower significance. (Without the presence of sky-lines, a derivation of the functional form of Ps{z) 
will be significantly more tractable with redshifts obtained from space-based spectroscopy, such as 
from Swift.) The large-scale rate results are much more sensitive to Pl{z), which is unfortunately 
the most difficult of the relevant probabilities to determine ab initio; in general Pl{z) should be 
constructed on a case-by-case basis, including any biases of distance upon triggering, localization 
probabilities, and redshift determination. There are certain psychological elements in the redshift 
determination (e.g., integrating longer until a redshift is found) that make this component partic- 
ularly difficult to model. It is clear, since the large-scale KS probability drops precipitously for 
zi ^ 1.25, that there must be some strong diminution of the detection rate of redshifts at larger 
distances. 

Given these difficulties, since we cannot test which star-formation rate GRBs trace best, we 
use the notion that GRBs should trace the universal SFR as a point of departure. That is, we 
fix a functional form of P{z) that provides a reasonable agreement with the expectations for the 
GRB rate on a large scale. Under this assumption (L = —2 and zi = 1), we test the significance 
of the apparent small-scale clustering. As can be seen in Table 3, there does not appear to be any 
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significant small-scale clustering of GRBs in redshift space when compared with the Monte Carlo 
set; however, the observed number of bursts paired with lA^I ^ 1000 km occurs in only one 
out of five iterations. 

The small-scale comparison in Table 3 is rather conservative, however. The only obvious a 
posteriori injection to the comparison is in choosing the velocity bins based upon the observed 
dataset. We have not required close pairs in the comparison set to fall within a certain redshift 
range, say z = [0.65, 2.05], nor have we required that the close pairs in the comparison sets be spaced 
by at least Sz = 0.1 (as observed). Any such restrictions would tend to increase the significance of 
the observed small-scale clustering. However, without some a priori hypothesis that the particular 
redshifts and spacings in the observed datasets are of interest, further restrictions are unwarranted. 

In the context of cosmologies with oscillating Hubble "constants" (see §1) some a priori pref- 
erence for special redshifts may indeed exist. For example, wc might expect that the most massive 
clusters of galaxies would reside at redshifts where the Hubble constant loiters at a local mini- 
mum. One of the most, if not the most massive clusters known, with the highest observed X-ray 
temperature and high velocity dispersion of member galaxies (cr > 1100 km s^^), is MS 1054—03 
(Neumann &; Arnaud 2000; van Dokkum et al. 2000). Correcting the systemic redshift reported in 
van Dokkum et al. (2000) to the local Hubble frame (following §2.4), the redshift of the cluster is 
zihf = 0.8337 ±0.0007. This is just 177 ± 125 km s"^ from the redshift of GRB 970508 yet the clus- 
ter and the GRB site are separated by = 88.4° on the sky. Further, the redshift of MS 1054—03 
is 1456 and 1872 km s"^ from GRB 990705 and GRB 000210 (respectively), less than two times 
the velocity dispersion of the cluster itself. Using our Monte Carlo simulation with Pg{z) oc SFl, 
the probability that three bursts out of 26 would fall within 2000 km of z = 0.8337 by random 
chance is 0.035. Since the association with a massive cluster was chosen a posteriori we do not 
claim this to be a significant result, but it should be of interest to test the association of new GRB 
redshifts with other massive clusters at moderate to high redshifts. 

The existence and expectation of several close groupings in redshift space holds some practical 
implications for observations of future GRBs. First, precise redshift determinations via moderate- 
resolution optical spectroscopy {R ^ 1000) will continue to be important even though approximate 
redshifts, using GRB distance-indicators (such as the Lag-Luminosity or Variability relations), 
may someday be used reliably to constrain the overall large-scale redshift distribution. Second, 
new future close redshift groupings will enable efficient detailed narrow-band imaging studies of 
multiple GRB hosts, such as the study of GRB 000926 and GRB000301C (Fynbo et al. 2002). 

We end by noting the curious absence of detected GRB redshifts in the redshift range between 
z ~ 2.3 — 3.2, where a relatively clean window for Lyman a emission detection exists at optical 
wavelengths (see fig. 1). One explanation for this is that redshift discovery via Ly a emission should 
be fruitless for about half of high redshift GRBs, since the average equivalent width of Ly a for 
Lyman Break Galaxies beyond 2; ~ 3 is near zero (K. Adelberger, personal communication); the 
absence of bursts in this range could simply be due to the low number of GRBs foUowed-up with 
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sufficiently small delay to detect Lyman a in absorption at optical wavelengths. In this respect, 
the uniformity and sheer rate of redshift discovery with Swift, should be most enlightening. 
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Fig. 1. — The distribution of recession velocity differences between nearest redshift pairs, following 
Table 2. The redshifts have been correction for the heliocentric motion through the local Hubble 
frame. Duari et al. 1992 performed a similar correction for quasar redshifts but only accounted for 
the heliocentric motion in the Galaxy. Representative velocity dispersions of galaxies and clusters 
are indicated with the dashed and dotted lines. There are 4 pairs of bursts which are separated by 
^ 1000 km s-^ 
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Fig. 2. — Comparison of the large-scale redshift distribution (solid cumulative line) with various models for 
redshift detection. Models (shown as dotted lines) that do not correct any observational biases in measuring 
high redshifts {L = 0) are clearly ruled out, but very different models for the true rate, when including a 
simple model for the high redshift bias (L = —2, zi ^ 1-25), are all allowed by the data. The Kolmogorov- 
Smirnov (KS) probabilities are given for six model comparisons, with the redshift of maximum deviation 
from the model noted with a vertical line. Clearly, the universal SFR cannot be inferred from the current 
sample. Here we have used GRB rate models that follow different parameterizations of the SFR (following 
Porciani & Madau 2001). SFl is the so-called Madau rate, where the GRB rate falls beyond redshift of 
unity. SF2 is a dust-corrected form of the Madau rate that levels off beyond redshift of unity (e.g., Steidel 
et al. 1999). SF3 is a rate which continues to increase beyond z !v 1.5. 
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Fig. 3. — A toy model for the redshift observabihty with ground-based optical spectroscopy in the 
presence of sky emission lines, Ps{z), for an instrumental resolution of 5 A. The observability of 
prominent star- formation emission lines (e.g., Ly a, [O II]) and gaseous halo absorption lines (e.g., 
Fe I, Mg II) have been modeled using the wavelengths of bright night sky lines as a mask; see 
text. The decrement of probability between 1.5 < z < 2 is due to the inaccessibility of redshiftcd 
Lya, Ha, and [O II] AA 3727 A in the optical bandpass. No optical redshift can be determined for 
bursts beyond z ^ 7, due to Lyman a blanketing of the optical spectra. The assignment of relative 
probability goes as follows: If one prominent emission line is observable, then we set ^5(2:0) = 0-9- 
If more than 1 emission line is observable then Ps{zo) = 1-0. If no emission lines are observable but 
at least a few prominent metal absorption lines are observable, then Ps{zq) = 0.3. If four of these 
nine absorption lines are observable then Ps{zq) = 0.4; for five observable lines, Ps{zq) = 0.5, etc. 
The main results of this paper are not strongly sensitive to the precise values of these (arbitrary) 
probability assignments. 
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Fig. 4. — A model distribution, P{z), of the relative probability for GRB redshift discovery using 
ground-based optical spectroscopy. The model is shown with low resolution for clarity. The overall 
shape is determined by the rate density, Pg{z), constructed assuming that the GRB rate should 
follow the Madau star-formation rate (SFl). Beyond redshift of z = 1 we assume that the relative 
detectability of emission features in GRB hosts, drop as the inverse square of the luminosity distance 
(L = —2). The uncertainty in the precise value of L (and hence Pl{z)) is also confounded by an 
uncertainty in the correct form of the rate density. The toy model for the redshift observability 
due to sky-lines is shown in Figure 3. The observed redshifts are noted as vertical lines above 
the distribution. A machine-readable data file of the high resolution version of this model may be 
obtained from the author by request. 
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Fig. 5. — Thresholds for the detection of significant small-scale clustering as a function of GRB 
redshift sample size. Based upon the Monte Carlo simulation (with SFl), the curves show the 
observed percentage of bursts required to be close to another burst as a function of velocity difference 
for 95% and 99% confidence. For instance in a sample of 100 ground-based optical redshifts, 66 
bursts must lie within 1000 km/s of another burst (i.e., 33 close pairs) to be considered a statistically 
significant excess. Note that for all samples sizes the percent requirements are similar for small 
velocity offsets, but diverge at large velocity offsets: with more bursts, the redshift density increases, 
so random pairing becomes more likely. 
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Table 1. Measured and corrected GRB redshifts 
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22 




+18:55:41 












971214 


11:56:26.0 


72.6 


109.1 


cm 


Lyct 


3.42(1) 


3.42(1) 


23 




+65:12:00 












000131 


06:13:31.0 


75.2 


93.2 


abs 


Lyman a forest 


4.511(2) 


4.513(2) 


24 




-51:56:40 















Note. — Column (5) gives the method of redshift determination for the burst (em = emission-line, abs = absorption-line) 
and column (6) lists the specific atomic species used to measure the redshift given in column (7). The redshift, corrected to the 
local Hubble frame, is given in column (8), following from equation 1. 

References. — 1. Garnavich ct al. (2003); 2. Vrccswijk ct al. (2001); 3. Price ct al. (2002b); 4. Price ct al. (2002c); 5. Bloom 
et al. (2001); 6. Djorgovski ct al. (1999); 7. Bloom ct al. (1998); 8. Lc Floc'h ct al. (2002); 9. Piro ct al. (2002); 10. Djorgovski 
et al. (2001a); 11. Djorgovski ct al. (1998); 12. Vrccswijk ct al. (1999); 13. Price ct al. (2002a); 14. Djorgovski ct al. (2000); 
15. Bloom et al. (2003); 16. Barth et al. (2003); 17. Castro et al. (2001); 18. Kulkarni et al. (1999); 19. Castro et al. (2000a); 
20. Castro et al. (2000b); 21. Holland et al. (2002); 22. IVIatheson et al. (2002); 23. Kulkarni et al. (1998); 24. Andersen et al. 
(2000) 



Table 2. Nearest GRB neighbors in redshift and velocity space 



Burst 


•Zlhf 


1 'I Cdl CO t 


Aft 


lAzI 


\Aii\ 


Name 




JDUlol IN dill tJ 


o 




km s 


GRB 


011121 


0.363(1) 


GRB 990712 


30.3 


0.070(1) 


15053 ± 223 


GRB 


990712 


0.4328(2) 


GRB 010921 


114.4 


0.0166(4) 


3457 ± 93 


GRB 


010921 


0.4494(4) 


GRB 990712 


114.4 


0.0166(4) 


3457 ± 93 


GRB 


020405 


0.69130(8) 


GRB 970228 


133.4 


0.0036(3) 


628 ± 55 


GRB 


970228 


0.6949(3) 


GRB 020405 


133.4 


0.0036(3) 


628 ± 55 


GRB 


991208 


0.7056(5) 


GRB 970228 


121.4 


0.0107(5) 


1887 ± 102 


GRB 


970508 


0.8348(3) 


GRB 990705 


152.1 


0.0078(3) 


1278 ± 59 


GRB 


990705 


0.8426(2) 


GRB 000210 


39.0 


0.0026(2) 


416 ± 47 


GRB 


000210 


0.8452(2) 


GRB 990705 


39.0 


0.0026(2) 


416 ± 47 


GRB 


970828 


0.9573(1) 


GRB 980703 


81.4 


0.0066(2) 


1008 ± 38 


GRB 


980703 


0.9639(2) 


GRB 970828 


81.4 


0.0066(2) 


1008 ± 38 


GRB 


991216 


1.02(1) 


GRB 000911 


42.3 


0.04(1) 


5419 ± 2967 


GRB 


000911 


1.0568(1) 


GRB 991216 


42.3 


0.04(1) 


5419 ± 2967 


GRB 


980613 


1.0974(2) 


GRB 000418 


54.6 


0.0229(2) 


3253 ± 35 


GRB 


000418 


1.1203(1) 


GRB 980613 


54.6 


0.0229(2) 


3253 ± 35 


GRB 


020813 


1.253(2) 


GRB 990506 


104.1 


0.057(2) 


7446 ± 266 


GRB 


990506 


1.309180(135) 


GRB 020813 


104.1 


0.057(2) 


7446 ± 266 


GRB 


010222 


1.4778(2) 


GRB 990123 


6.2 


0.1233(8) 


14550 ± 95 


GRB 


990123 


1.6011(8) 


GRB 990510 


126.2 


0.0184(17) 


2109 ± 195 


GRB 


990510 


1.6195(15) 


GRB 990123 


126.2 


0.0184(17) 


2109 ± 195 


GRB 


000301C 


2.0340(3) 


GRB 000926 


23.8 


0.0026(7) 


259 ± 75 


GRB 


000926 


2.0366(7) 


GRB 000301C 


23.8 


0.0026(7) 


259 ± 75 


GRB 


011211 


2.144(1) 


GRB 000926 


105.4 


0.107(1) 


10383 ± 119 


GRB 


021004 


2.328(1) 


GRB 011211 


163.0 


0.185(1) 


17082 ± 132 


GRB 


971214 


3.42(1) 


GRB 000131 


134.1 


1.09(1) 


65197 ± 654 


GRB 


000131 


4.513(2) 


GRB 971214 


134.1 


1.09(1) 


65197 ± 654 
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Table 3. Testing the Significance of the Observed Small-scale Structure 



\Av\ # GRBs paired Prob. (Observed | Simulation, L = —2) Simple Prob. 
km s-i < |At;| Pg{z) oc SFl oc SF2 oc SF3 



259 


2 


0.530 


0.471 


0.438 


0.408 


416 


4 


0.303 


0.247 


0.215 


0.252 


628 


6 


0.206 


0.145 


0.120 


0.191 


1008 


8 


0.187 


0.126 


0.096 


0.171 


1278 


9 


0.192 


0.120 


0.088 


0.160 


1997 


10 


0.402 


0.275 


0.213 


0.110 


2109 


12 


0.225 


0.133 


0.093 


0.099 


3253 


14 


0.373 


0.222 


0.159 


0.034 


3457 


16 


0.199 


0.092 


0.060 


0.026 


5419 


18 


0.364 


0.199 


0.135 


0.002 


7446 


20 


0.401 


0.230 


0.165 


6e-5 



Note. — Columns (3)-(6) give the probability estimates of observing, by 
random chance, at least the observed number of GRBs (column 2) that are 
paired with another within a velocity of |Ai;| (column 1). Columns (3)-(5) give 
the results from the Monte Carlo simulations described in the text for different 
functional forms of the universal star-formation rate. The last column gives 
the simplistic probability calculation of observing exactly the number of pairs 
assuming uniform selection in velocity space (equation 4 with k = col. [2]/2). 
The approximation breaks down for large velocities because the probabilities 
of observing more than the given number of bursts becomes non-negligible and 
relativistic velocity subtraction is not taken in to account. 



